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Abstract. We analyze generalized CP symmetries of two-Higgs doublet models, extending 
them from the scalar to the fermion sector of the theory. We show that, with a single exception, 
those symmetries imply massless fermions. The single model which accommodates a fermionic 
mass spectrum compatible with experimental data possesses a remarkable feature. It displays a 
new type of spontaneous CP violation, which occurs not in the scalar sector responsible for the 
symmetry breaking mechanism but, rather, in the fermion sector. 

1 Introduction 

Although the Standard Model (SM) of electroweak interactions has had extraordinary exper- 
imental confirmation, the Higgs mechanism of the symmetry breaking responsible for giving 
mass to the particles remains largely untested. In particular, there is no fundamental reason 
why the SM should have only one Higgs. In fact, one of the simplest extensions of the SM is the 
Two-Higgs Doublet Model (THDM) proposed by Lee pQ, where there are five scalar particles 
(three neutral and one charged). One of the main features of this model, and the main reason 
Lee proposed it, is that it allows for the possibility of the vacuum of the theory spontaneously 
breaking the CP symmetry . 

The most general scalar potential in the THDM contains 14 real parameters, as opposed to 
the 2 in the SM potential. There is thus a great interest in reducing the number of free parameters 
in the THDM, through the imposition of some symmetry in the theory. One such class of 
symmetries are the so-called generalized CP (GCP) transformations - field transformations which 
combine CP with unitary transformations. As far as we know, GCP were first discussed in [2]. 
Their explicit use for quarks appeared in [3] and GCP in the scalar sector was initially developed 
by the Viena group in [4} 02 [6]. GCP symmetries are quite interesting in that they reduce 
immensely to number of free parameters in the potential. Recently [8] it was shown that, 
from the point of view of the scalar potential, these GCP symmetries fall unto three categories, 
designated in [8] by CP1, CP2 and CP3. The CP1 class corresponds to the standard CP 
transformation. The CP3 models, unlike CP1 and CP2, correspond to a continuous symmetry. 

In this paper we will show that CP2 cannot be extended from the scalar sector to the 
fermionic one without forcing at least one fermion to have zero mass. Likewise, it will be 
shown that, out of the infinite number of possible CP3 models, only one can be extended to the 



fermionic sector without inducing zero fermion masses. Thus we are able to drastically curtail 
the number of models which includes an acceptable fermion sector. 

Having found only one acceptable GCP THDM, we will then proceed to analyze its CP 
properties. We conclude that the model has a new type of spontaneous CP breaking. In this 
model, CP violation does not emerge from explicit CP breaking in the Yukawa sector, as it 
does in the SM; nor does it emerge from spontaneous CP violation in the scalar sector, induced 
by the vacuum of the model, as in the Lee model. Rather, in the model herein presented, the 
lagrangian is explicitly CP conserving and the scalar sector preserves CP, even after spontaneous 
symmetry breaking. Nonetheless, the vacuum expectation values generated by the scalar sector 
induce a CP violating phase in the quark sector. As far as we know, this is the first illustration 
in the literature of this peculiar type of spontaneous CP violation. 

This paper is structured as follows: in section [2] we present the scalar potential of the theory 
and review recent results on its possible symmetries, as well as the impact that GCP has on its 
parameters. In section[3]we show how GCP can be extended to the fermion sector, and we prove 
that, with a single exception, GCP always implies massless fermions; we also perform a fast fit 
to experimental constraints. In section 0] we discuss the CP properties of this model, and show 
that it has a new type of spontaneous CP breaking. We present our conclusions in section [5j 



2 The scalar potential 



Let us first review the scalar sector of the THDM, and what is known of its symmetries. In 
the THDM one introduces two hypercharge Y = 1/2 Higgs scalar doublets. The most general 
THDM scalar potential which is renormalizable and compatible with the gauge symmetries of 
the SM is given by 

V H = m?]*}*! + m2 2 $2$2 - m? 2 $|$2 + H.c. + iAi($|$i) 2 + ±A 2 (34$ 2 ) 2 
+ A 3 ($ t 1 $i)($|$ 2 ) + A 4 ($I$ 2 )($5$i) + [iA 5 ($I$ 2 ) 2 + A 6 ($I$i)($l$ 2 ) 



A 7 ($5$ 2 )($I$2) +H.c. 



and it involves 14 parameters. Here m^, and Ai,--- , A4 are real parameters, while, in 
general, m 2 2 , A5, A6 and A7 are complex. "H.c." stands for Hermitian conjugation. 

This large number of parameters may be reduced by requiring that the model be invariant 
for some type of symmetry. Namely, we may require invariance for field transformations of the 
form 



S a b&b, 



for an S matrix belonging to U(2), which are known as Higgs family symmetries. Another pos- 
sibility consists on the so-called generalized CP symmetries HI El E] , requiring invariance 
under the following field transformations: 



$ a -> X ab <f>* b , 



(3) 



with X £ U{2) as well. Each of these symmetries corresponds to a different model, with different 
phenomenology. One may wonder how many different choices of scalar potentials can one make 
in the THDM. Recently it has been shown that applying the symmetries with any possible 



1 Both field transformations discussed here obey a basic requirement, namely leaving the kinetic terms of the 
lagrangian unaffected. 



choices for S and X leads only to six classes of scalar potentials [8]. The impact of the six 
classes on the parameters are given in Table HJ for specific basis choices. This result is already 



Table 1: Impact of the symmetries on the coefficients of the Higgs potential in a specific basis. 



See Ref. [8] for more details. 
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extremely important: there are only six possible different types of physical models in the scalar 
sector. The first three may be obtained from family symmetries alone, and are not the subject 
of this work. The next three may be obtained from GCP symmetries alone, and we will study 
them in detail. Let us begin by explaining what the difference is between them. Since any basis 
of fields {$i,$2} has to give the same physical results, we will use that liberty of basis choice 
to simplify the GCP field transformations of Eq. ([3]). In fact it has been shown that there is 
always a basis of scalar fields, for which the most general GCP transformation matrix X may 
be brought to the form 

ce s e 
-so c e 

where < 9 < ir/2 [9]. Henceforth, c = cos, s = sin and the subindices indicate the angle. The 
classes CPl, CP2, and CP3, arise respectively from Xq, X^/2, and Xq for any one (or several) 
9 excluding and ir/2. CPl is the usual CP symmetry. CP2, like CPl, is a discrete symmetry, 
although it eliminates far more parameters in the potential than CPl does. CP3 corresponds to 
an infinite number of possibilities, since the angle 9 varies in a continuous interval. We will now 
analyze the extension of these symmetries to the fermion sector. Although we will only look at 
the quark sector, the same conclusions can be trivially extended to leptons. 



Xa 



(4) 



3 Yukawa interactions 

In the previous section we have seen that there are only three classes of models with GCP 
symmetries. However, that study considered only the scalar sector, not the whole lagrangian. 
We will now show how GCP can be extended to the fermion sector, and what consequences arise 
thereof. The scalar-quark Yukawa interactions of the THDM may be written as 

-C Y = q L (Ti^t + r 2 $ 2 )n H + ?l(Ai$i + A 2 ^ 2 )pr + H.c, (5) 

where qL = (pL, n L) T { n R an d pr) is a vector in the 3-dimensional generation space of left- 
handed doublets (right-handed charge —1/3 and +2/3) quarks. T\, T 2 , Ax, and A2 are 3x3 
matrices. 

We now wish to impose GCP on scalars and also on fermions. For the quark fields, the GCP 
transformations take the form 

n R ^X^°Cn* R , 

p R -+X 7 j°Cp R , (6) 



where 7 (C) is the Dirac (charge-conjugation) matrix. As was the case with the scalar GCP 
transformations, we can also simplify the above with an appropriate basis choice for the fermion 
fields. In fact, we may take X to the simple form [9] 



X, 














1 



(7) 



where < a < tt/2, and similarly for Xp and X-y. GCP is a good symmetry of Cy if and only 



ifr* 



(Xg) a bXaT a Xs, which we may write as 



X a Tl - (cgTi - s e T 2 )Xp = 0, 
X a T* 2 - (s e T l + c e T 2 )X p = 0. 



(8) 



Since they are complex, this gives us 36 equations in the 36 unknown real and imaginary parts 
of the various entries of the V matrices. Because of Eq. (J7]), they break into four blocks, which 
we designate by mn, m3, 3n, and 33, where m and n can take the values 1 or 2. In each block we 
have a system of homogeneous linear equations; the parameters are zero unless the determinant 
of the system vanishes. For instance, the real (Re) and imaginary (Im) parts of (1^1)33 and (1^)33 
in Eq. ([8]) give the following set of linear equations: 
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0, 
= 0. 



(9) 



The determinant of the first of these matrices gives 2(1 — eg), whereas the second one is equal to 
2(1 + eg). Because of the limited range of variation of 9, the second determinant is never zero - 
which means that the T33 coefficients will always be real. As for the first of these determinants, 
it is only zero if 9 = (the standard CP definition). For any other case, the 33 entries of the T 
matrices will be zero. 

A similar reasoning applies to the m2> entries. Eq. ([8]) gives us a set of 4 x 4 equations on the 
real and imaginary parts of the Yukawa couplings, with corresponding determinants given by 
4(c Q — eg) 2 and 4(c Q + eg) 2 , respectively. Since a has the same range of variation than 9, these 
determinants only vanish in very specific situations; they vanish when a = 9oia = 9 = tt/2, 
respectively. In the same manner, the 3n equations give rise to determinants equal to 4(cg — ce) 2 
and 4(cg + ce) 2 for the real and imaginary parts, respectively, with analogous conditions for 
vanishing determinants. 

Finally, the mn block from Eq. ([8]) gives us a set of 8 x 8 equations for the real and imaginary 
parts of the corresponding Yukawas. The corresponding determinants are 



16 (eg ± c a+p f (eg ± c Q 



(10) 



with the "+" signs corresponding to the equations for the imaginary parts. We performed a 
thorough analysis of the conditions for vanishing determinants, and the results are summarized 
in Table [2j Similar results hold for the charged +2/3 quark matrices Ai and A2, with /3 — >■ 7. 

After spontaneous symmetry breaking, the fields acquire the vacuum expectation values 
(vevs) vi/y/2 and V2e tS /y/2; where v\ and v 2 are real, without loss of generality. It is convenient 
to rotate into the so-called Higgs basis {Hi,H 2 } through <£ a = U a bHb, where 



1 

v 
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v 2 e 
—v\e~ 



(11) 



Table 2: Impact of the GCP symmetry on T\ and T2- We separate real (Re) and imaginary 
(Im) components. 
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(12) 



where G + and G° are the Goldstone bosons, which, in the unitary gauge, become the longitudinal 
components of the W + and of the Z°, and H°, R and / are real neutral fields. In the new scalar 
basis, the Yukawa coupling matrices become = TbUb a , A^ = AbU b * a , and the scalar couplings 
are also rotated. Finally, the fermion mass basis is obtained through transformations U a with 
a = di,dji,UL,UR that diagonalize the quark couplings to Hi, 



D d = diag(m d , m s ,m b ) 
% Ul L Af U UR = D u = diag(m M , m c ,m t ) 

while the couplings to Hi become (we follow the notation of |10|. [TT] ) 



(13) 



(v/V2) c/^rf U dR = N d , 
(v/V2)Ul L A?U UR = N u . 

The Yukawa lagrangian may then be written as 



(14) 



V2 



L Y = (u L V, d L )(D d Hi + N d H 2 ) d R 

+ (u L , d L V^)(D u Hi + N U H 2 ) u R + H.c. 



(15) 



where V = lll L U dL is the CKM matrix. Thus N d and N u are responsible for flavor changing 
neutral currents (FCNC) involving the scalars R and I. 

Let us now look back at Table [2j We see that 9 = 0, corresponding to the usual definition of 
CP, forces all Yukawa couplings to be real. This is well known, but it requires that CP violation 
arise spontaneously. This theory leads to FCNC, which are constrained by experiment. This 
could be cured with further discrete symmetries. It is known that imposing the absence of FCNC 
is inconsistent with GCP for three quark families and any number of Higgs fields [6]. Here we 



will take the view that the scalar masses may be large enough to suppress FCNC, and consider 
the most general couplings. 

The case of 9 = tt/2 corresponds to CP2. As we have seen, this forces (ri)33 = ^2)33 = 0. 
If a = /3 = 7r/2 then the mn block vanishes. Thus the determinant of D d vanishes, forcing one 
quark mass to zero. This is excluded by experiment. If a 7^ tt/2 and/or (3 7^ tt/2, then the last 
column and/or the last row of D d vanishes, which leads again to a zero mass quark. As a result, 
it is impossible to extend CP2 to the quark sector in a way consistent with experiment. This 
zero mass problem had already been faced in |12} 113] . for a particular extension of CP2. We 
note that CP2 could be consistently extended to the fermion sector provided there were four 
fermion families. 

Now we come to the most interesting case: consider < 9 < tt/2. Here (ri)33 = ^2)33 = 
and, in order for the m3 and 3n entries to differ from zero (otherwise there would be at least a 
zero mass quark), we need to have a = f3 = 9. In this case, the determinants for the mn block 
in Eq. (jlOp simplify to 



Imaginary: 256 (c e/2 ) 8 (l - 2c e f 
Real: 256 (s e/2 f(l + 2c e f . 



(16) 
(17) 



The second determinant never vanishes, which implies that the mn entries will only have imag- 
inary parts. And the first determinant, Eq. (|16p . only vanishes if 9 = ir/3. For all other cases, 
the quark mass matrices will have zero eigenvalues. Notice that the fact that the angles a, (3 
and 9 are constrained to be in the interval [0,7r/2] was fundamental in this demonstration. 

Therefore, out of the infinite number of GCP symmetries, only one survives: the case where 
one has 9 = a = (3 = ir /3. For this model the Yukawa matrices have the form: 



mn «ai2 tll3 

mi2 -mn ^23 
a?,i a 32 

iai2 -mn -023 
-ian -mi2 ai 3 
-032 031 



(18) 



Notice that Ti has 2 (4) independent imaginary (real) entries and that T2 does not involve 
any new parameter. Similar parametrizations hold for Ai and A2, involving 6 new parameters 
(a — > b). We stress that this is the only possible extension of GCP into the fermion sector 
consistent with the fact that the quarks have non-vanishing masses, and it leads to a very 
tightly constrained "minimal" model. 

Our theory has only 12 Yukawa parameters and two independent vevs. Indeed, v\ may be 
made real without loss of generality; given v 2 , this fixes v 2 , and we only need the phase 5. In 
order to perform a numerical fit, we define the matrices 



£lf (Tf )t = u dj 
i Af (Af )t = U UL DM 



(19) 



The eigenvalues of these matrices give the square of the quark masses, providing six constraints. 
There are four other experimental constraints arising from the four independent parameters 
parametrizing the CKM matrix V. We can take these as the eigenvalues of H u H d , and also the 
CP-violating quantity [3] 

Tr[tf n , H d f = 6i(m 2 t - m 2 c ){m 2 t - m 2 u )(m 2 c ' 



J 



m,. 



x (ml - m 2 )(m 2 - m\) (m 2 - m 2 d )lm (V us V cb V* b V* s ) 



(20) 



Using Eq. (fTBj) in Eq. (f20l) . we have shown that 5 = implies J = 0. Said otherwise, spontaneous 
CP violation is required in order to have CP violation in the CKM matrix, and this only occurs 
if we add a soft symmetry breaking term to the scalar potential. However, the CP properties of 
this model are far from obvious, and will be discussed in the following section. 

Using a fast numerical fit, we have found that setting an = 4.6927 x 1CU 6 , ai 2 = —5.9799 x 
10" 4 , ai3 = -2.32 x 1(T 2 , a 23 = -6.6 x 10~ 3 , a 31 = -8.815 x 10~ 5 , a 32 = 5.1193 x 1(T 6 
in the down-quark sector and b n = 7.3 x 1CT 3 , 6 12 = 7.6445 x 1(T 5 , b 13 = 9.578 x 10 -1 , 
623 = 2.325 x 10 -1 , 631 = 1.3446 x 10~ 4 , 6 32 = 5.9491 x 10" 4 , in the up-quark sector, and 
vi = 173.944, v 2 cosS = -0.8467, w 2 sin<5 = -0.9565, we obtain m d = 0.00298, m s = 0.10511, 
nib = 4.19701, m u = 0.00200, m c = 1.27439, mt = 171.451, for the masses, and the magnitudes 
of the CKM matrix elements are 



\V\ 



0.97430 0.22521 0.00339 
0.22516 0.97348 0.04039 
0.00579 0.04011 0.99918 



(21) 



All masses and vevs are in GeV. The values of v 2 = v 2 + v 2 , the masses, and the magnitudes 
of the CKM matrix elements are in very good agreement with the experimental data |14j . One 
may wonder whether the vevs we used above can be generated by the scalar potential. As we 
mentioned earlier, we need to introduce soft breaking terms in the scalar potential to generate 
5, namely by setting mf^ 7^ m| 2 and Re(mf 2 ) + 1 0. It is easy to find values for the potential's 
parameters which generate the desired vevs: they can be reproduced by m 2 ^ = —28791, m 22 = 
5679.5, Re(mf 2 ) = -167.8012, Im(m? 2 ) = 0, Ai = A 2 = 0.9516, A 3 = 0.0176, and A 4 = 0.3643 
(the quadratic parameters are expressed in units of GeV 2 ) . 

The numerical fit is quite successful, but the exception is \ Vtd\ which, due to the hierarchical 
nature of the CKM matrix elements, becomes much more difficult to fit H. A related problem 
occurs with the CP-violating quantity J c km = im (KisU c bU u * b U c * ). For our choice of parameters, 
one obtains J c km — — 5.9x 10 -8 , while the experimental data leads to J c km = (3.05+0.20) x 10 -5 . 
The reason has to do with the sensitivity of our fit procedure to the input parameters. We have 
fit Z us , Z u b, Z cs , and Z c b, where Zij = \Vij\ 2 . These four elements can be used to parametrize 
completely the CKM matrix [151 CE1 DZ] • In particular [18] , 

^ckm = ^ZugZubZcgZcb — (1 — Z us — Z u b — Z cs — Z c b + Z us Z c b + Z u hZ cs ) 2 . (22) 

Using the experimental values for the CKM matrix elements |14j . one can show that to first 
order 

^htm „ 540.95^ + 0.7^ 

Jckm Z us Z u b 

+ 10065.9—^ + 18.0—^. (23) 

Z C s Z c b 

This shows how extremely sensitive J c km is to the exact value of V cs , explaining why it is easy 
to perform a fast fit to the latter, but thus inducing a large error on the former. 



4 A New Type of Spontaneous CP Violation 

We now come to the most remarkable feature of this model, namely its spontaneous CP violation 
properties. We have already mentioned that, to obtain a non-zero value for the Jarlskog invariant 

2 Notice that the direct measurement of \Vtd\ comes from the mixing in the B system, which occurs in the SM 
through a box diagram. This mixing can receive contributions involving the second Higgs, thus altering the value 
of IVwI. 



J, Eq. ()20p . we need the vevs of the scalar fields to have a relative phase. Usually, this fact 
by itself is taken to be a sign of spontaneous CP violation in the scalar sector. The model 
we are now discussing behaves, however, in a very different fashion. Regardless of its possible 
future experimental relevance as a viable description of nature, this model serves as a theoretical 
example of a new kind of spontaneous CP violation. 

First, let us demonstrate that the lagrangian of this model is explicitly CP-conserving, i.e. 
that it is possible to write it in a field basis where all parameters are real. We see in the CP3 
entry in Table Q] that the scalar potential only has real parameters. As was mentioned above, 
the soft breaking terms we introduced are also real, so there are no complex phases in the scalar 
sector, prior to SESB. 

Now for the Yukawa sector. It might seem, looking at the form of the Yukawa matrices in 
Eq. (|18p . that one cannot avoid the presence of complex numbers in the Lagrangian. However, 
we are free to choose a different quark basis, by performing transformations of the type 



QL -> U L q L , n R -> UrUr, 
with ?7(3) matrices Ul and Ur. This changes the Yukawa matrices according to 

r k -> ulr k u R . 



(24) 



(25) 



By choosing lit = Ur = diag(e l?r//4 , e e l7f//4 ), we remove all factors of i from the V matrices 



and are left with 



an 
ai2 
031 

a 12 
-an 
-CJ32 



ai2 «13 
-an 023 
032 



-an 
-012 

«31 



-023 

013 




(26) 



where all parameters are real. We stress that the basis choice that leads to this result involves 
only the fermion fields, and as such does not introduce any phases in the scalar potential. An 
identical fermion basis choice may be used to render real the Yukawa matrices in the up-quark 
and leptonic sectors. Hence, we have proved that it is possible to find a basis where all parameters 
in the lagrangian are real - the model is explicitly CP conserving. 

After SESB, with the soft breaking we discussed, we introduce a complex phase 5 in the 
theory, from the vevs Vi/y/2 and V2e lS /V2. One might wonder whether 5^0 implies the 
presence of spontaneous CP violation in the scalar sector. This is best investigated with the 
basis-invariant quantities developed by Lavoura and Silva [19] and by Botella and Silva 
Their explicit calculation shows that they all vanish - thus, with the proposed soft breaking of 
the CP3 symmetry there is no CP violation in the scalar sector, whether explicit or spontaneous. 
A simple way to confirm this is by changing to the Higgs basis, as given by the transformation 
in Eq. (jlip . In that basis, the potential, with the soft breaking terms already included, becomes 



m n #i#i + rh\ 2 n\ll 2 



in 



12 J 



H\H 2 + H.c. + fX^HlH^ 2 + f\ 2 (HlH 2 ) 



+ \z(HlH x )(HlH 2 ) + \ A (HlH 2 )(HlH x ) + f\ 5 (HlH 2 ) 2 + X^HlH^Hl^ 



+ A 7 (i/ 2 t ^ 2 )(^i t ^2)+H.c. 



(27) 



O 




where 




) 



(29) 



The only phases invariant under a trivial rephasing of H2 and, thus, possibly signaling CP 
violation are In^AgA?), Im(A|Ag), and Im(A| Ay) . These are precisely the three quantities 
introduced by Lavoura and Silva [19J to probe all possible sources of CP violation in the scalar 
sector in a basis-invariant way. Since these quantities are all zero, we conclude that there is no 
CP violation in the scalar sector of the theory, even after SESB. This means, in particular, that 
one can find a field basis where there is no mixing between the CP-even and CP-odd neutral 
scalar particles. 

One may then ask whether the phase 5 is at all relevant in the model. The answer is yes 
because, due to the interplay between the scalar and Yukawa sectors of the theory, there is no 
choice of field basis through which one can absorb 5. In fact, the CP-violating quantity J of 
Eq. (|20l) is directly proportional to sin 5. Notice, too, that as long as 5 ^ the FCNC involving 
Nd and N u , in Eq. (|14p . will also, in general, involve complex phases and as such may well serve 
as further sources of CP violation^ But all possible sources of CP violation vanish if 5 = 0, 
even though no CP breaking occurs in the scalar sector. 

This, then, is a new type of CP violation: 

• It is not like that of the SM, since there CP is explicitly broken at the lagrangian level; 

• It is not like the CP violation that occurs in Lee-type mechanisms, since there CP is 
spontaneously broken in the scalar sector. 

In the model herein presented, it is the fermion sector which exhibits the CP violation that 
arises spontaneously. And the CKM matrix is generated through a spontaneous breaking of CP, 
not an explicit one as is the case of the SM. Of course, these results were obtained at tree level. 
Still, it is interesting that the scalar sector does the deed (spontaneously break the symmetry) 
but it is the fermion sector which pays the consequence (providing CP violation). To the best 
of our knowledge, this type of CP violation is unheard of in the literature. 

5 Conclusions 

We have shown that, for the most part, the generalized CP symmetries of the scalar sector cannot 
be extended to the quark sector with three generations, while keeping all six quarks massive. 
This reduces dramatically the types of available models. We have found the only exception: the 
only THDM with GCP which leads to an acceptable fermion mass spectrum. We have shown 
that the model has very few parameters, both in the scalar and Yukawa sectors, and that it 
can give a good fit to the known quark masses and mixings. Due to the hierarchical nature 

3 If the masses of the new scalars are large enough, then these sources of CP violation will have a small impact 
on current experiments, including those currently used to constrain V . Our aim in this section is to highlight a 
new scenario for spontaneous CP violation and not a particular implementation thereof. 



n 



of the masses and CKM matrix elements, very precise numerical fits are needed. This model 
has FCNC which might have considerable impact on its experimental predictions for hadron 
phenomenology. A detailed study of its impact on such observables is therefore needed, and will 
be presented elsewhere |20| . 

We have also shown that this model has the peculiar feature that the scalar sector by itself 
is both explicitly and spontaneously CP-conserving, but that it generates a symmetry breaking 
which induces CP violation in the quark sector. This constitutes a new type of spontaneous 
CP violation, different from the SM or the usual findings of the THDM. As such this model, 
regardless of its experimental relevance, serves as the first example of a new way of generating 
CP violation. It is not unique in that respect: in fact, it is a trivial exercise to build a three-Higgs 
doublet model which displays exactly the same CP properties. 

There is a further motivation for studying in detail models of spontaneous CP violation. 
The excellent agreement between the SM and all CP violating experiments has been taken as 
a definitive confirmation of the SM's CP violation mechanism, i.e. the existence of explicitly 
CP violating Yukawa couplings. Models of CP violation, like the one presented here, where J 
could be the dominant source of CP violation raise the tantalizing possibility that, contrary 
to widespread belief, the current experiments on CP violation confirm the Cabibbo-Kobayashi- 
Maskawa source of CP violation but not its origin in the explicit CP breaking of the Yukawa 
interactions. 
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